A danger of low copy numbers for inferring incorrect cooperativity degree

Background: A dose-response curve depicts fraction of bound proteins as a function of unbound ligands. Dose-response curves are used to measure the cooperativity degree of a ligand binding process. Frequently, the Hill function is used to fit the experimental data. The Hill function is parameterized by the value of the dissociation constant, and the Hill coefficient which describes the cooperativity degree. The use of Hill's model and the Hill function have been heavily criticised in this context, predominantly the assumption that all ligands bind at once, which lead to further refinements of the model. In this work, the validity of the Hill function has been studied from an entirely different point of view. In the limit of low copy numbers the dynamics of the system becomes noisy. The goal was to asses the validity of the Hill function in this limit, and to see in which ways the effects of the fluctuations change the form of the dose-response curves. Results: Dose-response curves were computed taking into account effects of fluctuations. The effects of fluctuations were described at the lowest order (the second moment of the particle number distribution) by using previously developed Pair Approach Reaction Noise EStimator (PARNES) method. The stationary state of the system is described by nine equations with nine unknowns. To obtain fluctuation corrected dose-response curves the equations have been investigated numerically. Conclusions: The Hill function cannot describe dose-response curves in a low particle limit. First, dose-response curves are not solely parameterized by the dissociation constant and the Hill coefficient. In general, the shape of a dose-response curve depends on the variables that describe how an experiment (ensemble) is designed. Second, dose-response curves are multi valued in a rather non-trivial way

A new multicompartmental reaction-diffusion modeling method links transient membrane attachment of E. coli MinE to E-ring formation

Many important cellular processes are regulated by reaction-diffusion (RD) of molecules that takes place both in the cytoplasm and on the membrane. To model and analyze such multicompartmental processes, we developed a lattice-based Monte Carlo method, Spatiocyte that supports RD in volume and surface compartments at single molecule resolution. Stochasticity in RD and the excluded volume effect brought by intracellular molecular crowding, both of which can significantly affect RD and thus, cellular processes, are also supported. We verified the method by comparing simulation results of diffusion, irreversible and reversible reactions with the predicted analytical and best available numerical solutions. Moreover, to directly compare the localization patterns of molecules in fluorescence microscopy images with simulation, we devised a visualization method that mimics the microphotography process by showing the trajectory of simulated molecules averaged according to the camera exposure time. In the rod-shaped bacterium _Escherichia coli_, the division site is suppressed at the cell poles by periodic pole-to-pole oscillations of the Min proteins (MinC, MinD and MinE) arising from carefully orchestrated RD in both cytoplasm and membrane compartments. Using Spatiocyte we could model and reproduce the _in vivo_ MinDE localization dynamics by accounting for the established properties of MinE. Our results suggest that the MinE ring, which is essential in preventing polar septation, is largely composed of MinE that is transiently attached to the membrane independently after recruited by MinD. Overall, Spatiocyte allows simulation and visualization of complex spatial and reaction-diffusion mediated cellular processes in volumes and surfaces. As we showed, it can potentially provide mechanistic insights otherwise difficult to obtain experimentally

Graphical Approach to Model Reduction for Nonlinear Biochemical Networks

Model reduction is a central challenge to the development and analysis of multiscale physiology models. Advances in model reduction are needed not only for computational feasibility but also for obtaining conceptual insights from complex systems. Here, we introduce an intuitive graphical approach to model reduction based on phase plane analysis. Timescale separation is identified by the degree of hysteresis observed in phase-loops, which guides a “concentration-clamp” procedure for estimating explicit algebraic relationships between species equilibrating on fast timescales. The primary advantages of this approach over Jacobian-based timescale decomposition are that: 1) it incorporates nonlinear system dynamics, and 2) it can be easily visualized, even directly from experimental data. We tested this graphical model reduction approach using a 25-variable model of cardiac β1-adrenergic signaling, obtaining 6- and 4-variable reduced models that retain good predictive capabilities even in response to new perturbations. These 6 signaling species appear to be optimal “kinetic biomarkers” of the overall β1-adrenergic pathway. The 6-variable reduced model is well suited for integration into multiscale models of heart function, and more generally, this graphical model reduction approach is readily applicable to a variety of other complex biological systems

Temperature Control of Fimbriation Circuit Switch in Uropathogenic Escherichia coli: Quantitative Analysis via Automated Model Abstraction

Uropathogenic Escherichia coli (UPEC) represent the predominant cause of urinary tract infections (UTIs). A key UPEC molecular virulence mechanism is type 1 fimbriae, whose expression is controlled by the orientation of an invertible chromosomal DNA element—the fim switch. Temperature has been shown to act as a major regulator of fim switching behavior and is overall an important indicator as well as functional feature of many urologic diseases, including UPEC host-pathogen interaction dynamics. Given this panoptic physiological role of temperature during UTI progression and notable empirical challenges to its direct in vivo studies, in silico modeling of corresponding biochemical and biophysical mechanisms essential to UPEC pathogenicity may significantly aid our understanding of the underlying disease processes. However, rigorous computational analysis of biological systems, such as fim switch temperature control circuit, has hereto presented a notoriously demanding problem due to both the substantial complexity of the gene regulatory networks involved as well as their often characteristically discrete and stochastic dynamics. To address these issues, we have developed an approach that enables automated multiscale abstraction of biological system descriptions based on reaction kinetics. Implemented as a computational tool, this method has allowed us to efficiently analyze the modular organization and behavior of the E. coli fimbriation switch circuit at different temperature settings, thus facilitating new insights into this mode of UPEC molecular virulence regulation. In particular, our results suggest that, with respect to its role in shutting down fimbriae expression, the primary function of FimB recombinase may be to effect a controlled down-regulation (rather than increase) of the ON-to-OFF fim switching rate via temperature-dependent suppression of competing dynamics mediated by recombinase FimE. Our computational analysis further implies that this down-regulation mechanism could be particularly significant inside the host environment, thus potentially contributing further understanding toward the development of novel therapeutic approaches to UPEC-caused UTIs

Leveraging analytics to produce compelling and profitable film content

Producing compelling film content profitably is a top priority to the long-term prosperity of the film industry. Advances in digital technologies, increasing availabilities of granular big data, rapid diffusion of analytic techniques, and intensified competition from user generated content and original content produced by Subscription Video on Demand (SVOD) platforms have created unparalleled needs and opportunities for film producers to leverage analytics in content production. Built upon the theories of value creation and film production, this article proposes a conceptual framework of key analytic techniques that film producers may engage throughout the production process, such as script analytics, talent analytics, and audience analytics. The article further synthesizes the state-of-the-art research on and applications of these analytics, discuss the prospect of leveraging analytics in film production, and suggest fruitful avenues for future research with important managerial implications

Improved functionalization of oleic acid-coated iron oxide nanoparticles for biomedical applications

Superparamagnetic iron oxide nanoparticles can providemultiple benefits for biomedical applications in aqueous environments such asmagnetic separation or magnetic resonance imaging. To increase the colloidal stability and allow subsequent reactions, the introduction of hydrophilic functional groups onto the particles’ surface is essential. During this process, the original coating is exchanged by preferably covalently bonded ligands such as trialkoxysilanes. The duration of the silane exchange reaction, which commonly takes more than 24 h, is an important drawback for this approach. In this paper, we present a novel method, which introduces ultrasonication as an energy source to dramatically accelerate this process, resulting in high-quality waterdispersible nanoparticles around 10 nmin size. To prove the generic character, different functional groups were introduced on the surface including polyethylene glycol chains, carboxylic acid, amine, and thiol groups. Their colloidal stability in various aqueous buffer solutions as well as human plasma and serum was investigated to allow implementation in biomedical and sensing applications.status: publishe

Search for Dark Matter and Supersymmetry with a Compressed Mass Spectrum in the Vector Boson Fusion Topology in Proton-Proton Collisions at root s=8 TeV

Peer reviewe

Search for strongly interacting massive particles generating trackless jets in proton-proton collisions at s = 13 TeV

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb - 1 , collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses

Search for dark matter produced in association with a single top quark or a top quark pair in proton-proton collisions at s=13 TeV

A search has been performed for heavy resonances decaying to ZZ or ZW in 2l2q final states, with two charged leptons (l = e, mu) produced by the decay of a Z boson, and two quarks produced by the decay of a W or Z boson. The analysis is sensitive to resonances with masses in the range from 400 to 4500 GeV. Two categories are defined based on the merged or resolved reconstruction of the hadronically decaying vector boson, optimized for high- and low-mass resonances, respectively. The search is based on data collected during 2016 by the CMS experiment at the LHC in proton-proton collisions with a center-of-mass energy of root s = 13 TeV, corresponding to an integrated luminosity of 35.9 fb(-1). No excess is observed in the data above the standard model background expectation. Upper limits on the production cross section of heavy, narrow spin-1 and spin-2 resonances are derived as a function of the resonance mass, and exclusion limits on the production of W' bosons and bulk graviton particles are calculated in the framework of the heavy vector triplet model and warped extra dimensions, respectively.A search has been performed for heavy resonances decaying to ZZ or ZW in 2l2q final states, with two charged leptons (l = e, mu) produced by the decay of a Z boson, and two quarks produced by the decay of a W or Z boson. The analysis is sensitive to resonances with masses in the range from 400 to 4500 GeV. Two categories are defined based on the merged or resolved reconstruction of the hadronically decaying vector boson, optimized for high- and low-mass resonances, respectively. The search is based on data collected during 2016 by the CMS experiment at the LHC in proton-proton collisions with a center-of-mass energy of root s = 13 TeV, corresponding to an integrated luminosity of 35.9 fb(-1). No excess is observed in the data above the standard model background expectation. Upper limits on the production cross section of heavy, narrow spin-1 and spin-2 resonances are derived as a function of the resonance mass, and exclusion limits on the production of W' bosons and bulk graviton particles are calculated in the framework of the heavy vector triplet model and warped extra dimensions, respectively.A search for dark matter produced in association with top quarks in proton-proton collisions at a center-of-mass energy of 13 TeV is presented. The data set used corresponds to an integrated luminosity of 35.9 fb(-1) recorded with the CMS detector at the LHC. Whereas previous searches for neutral scalar or pseudoscalar mediators considered dark matter production in association with a top quark pair only, this analysis also includes production modes with a single top quark. The results are derived from the combination of multiple selection categories that are defined to target either the single top quark or the top quark pair signature. No significant deviations with respect to the standard model predictions are observed. The results are interpreted in the context of a simplified model in which a scalar or pseudoscalar mediator particle couples to a top quark and subsequently decays into dark matter particles. Scalar and pseudoscalar mediator particles with masses below 290 and 300 GeV, respectively, are excluded at 95% confidence level, assuming a dark matter particle mass of 1 GeV and mediator couplings to fermions and dark matter particles equal to unity.Peer reviewe

Search for the pair production of light top squarks in the e(+/-)mu(-/+) final state in proton-proton collisions at root s=13 TeV

A search for the production of a pair of top squarks at the LHC is presented. This search targets a region of parameter space where the kinematics of top squark pair production and top quark pair production are very similar, because of the mass difference between the top squark and the neutralino being close to the top quark mass. The search is performed with 35.9 fb(-1) of proton-proton collisions at a centre-of-mass energy of root s = 13 TeV, collected by the CMS detector in 2016, using events containing one electron-muon pair with opposite charge. The search is based on a precise estimate of the top quark pair background, and the use of the M-T2 variable, which combines the transverse mass of each lepton and the missing transverse momentum. No excess of events is found over the standard model predictions. Exclusion limits are placed at 95% confidence level on the production of top squarks up to masses of 208 GeV for models with a mass difference between the top squark and the lightest neutralino close to that of the top quark.Peer reviewe